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Abstract. Prudent walks are special self-avoiding walks that never take a step towards an already occupied 
site, and k-sided prudent walks (with k = 1,2,3,4) are, in essence, only allowed to grow along k directions. 
Prudent polygons are prudent walks that return to a point adjacent to their starting point. Prudent walks 
and polygons have been previously enumerated by length and perimeter (Bousquet-Melou, Schwerdtfeger; 
2010). We consider the enumeration of prudent polygons by area. For the 3-sided variety, we find that the 
generating function is expressed in terms of a q-hypergeometric function, with an accumulation of poles 
towards the dominant singularity. This expression reveals an unusual asymptotic structure of the number of 
polygons of area n, where the critical exponent is the transcendental number log 2 3 and and the amplitude 
involves tiny oscillations. Based on numerical data, we also expect similar phenomena to occur for 4-sided 
polygons. The asymptotic methodology involves an original combination of Mellin transform techniques and 
singularity analysis, which is of potential interest in a number of other asymptotic enumeration problems. 



1. Introduction 

The problem of enumerating self-avoiding walks (SAWs) and polygons (SAPs) on a lattice is a famous 
one, whose complete solution has thus far remained most elusive. For the square lattice, it is conjectured 
that the number SAW„ of walks of length n and the number SAP,, of polygons of perimeter n each satisfy 
an asymptotic formula of the general form 

(1) C-pT-nP, 

where C, /.i € M>n and (3 € R. (In the case of polygons, it is understood that n must be restricted to even 
values.) The number fj, is the "growth constant" and the number (3 is often referred to as the "critical 
exponent" . More precisely, the following expansions are conjectured, 

(2) SAW„ ~ Ci • fi n - n 11 / 32 , SAP„ ~ C 2 ■ //"■ n~ 5/2 , 

n— >oo 7i— >oo 

for some Ci, C2 > 0. 

For the square lattice, numerical methods based on acceleration of convergence and differential approx- 
imants suggest the value /1 = 2.6381585303 .... This estimate is indistinguishable from the solution of the 
biquadratic equation 13/x 4 — 7/i 2 — 581 = 0, which we consider to be a useful mnemonic. This was observed by 
Conway, Enting and Guttmann [5] in 1993, and verified to 11 significant digits by Jensen and Guttmann |23j 
in 2000, based on extensive numerical analysis of the sequence (SAP2„) up to In = 90. (Remarkably enough, 
for the honeycomb lattice, it had long been conjectured that the growth constant of walks is the biquadratic 
number a fact rigorously established only recently by Duminil-Copin and Smirnov fT3"j.) 

As regards critical exponents, the conjectured value (3 = |^ for walks is supported by results of Lawler, 
Schramm, and Werner that relate the self-avoiding walk to the "stochastic Loewner Evolution" (SLE) process 
of index 8/3; see, for instance, the account in Werner's inspiring lecture notes |37j . For (unrooted) polygons, 
the value /3 = — | was suggested by numerical analysis of the exact counting sequence, with an agreement 
to the seventh decimal place [53]. It is also supported by the observation that many simplified, exactly 
solvable, (naturally rooted) models of self-avoiding polygons appear to exhibit an n~ 3 / 2 universal behaviour 
- for these aspects, we refer to the survey by Bousquet-Melou and Brak 0], as well as the books [TBI |2"UI l3"o] . 

Regarding lattice polygons, which are closed walks, there is also interest in enumeration according to 
area, rather than perimeter. This question has analogies with the classical unsolved problem of enumerating 
polyominoes, also known as animals (these may have "holes"), according to the number of cells they contain. 
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Conjccturally [24) . the number a n (respectively, b n ) of polygons (respectively, polyominoes) comprised of n 
cells satisfy asymptotic estimates of the form 

a n ~ C 3 • (3.9709...)™ -n' 1 , b n ~ C 4 • (4.0625 ...)"• rT 1 , 

n— j-cjo n— >oo 

for some 6*3,(74 € K>o; these asymptotic estimates are still of the form (JTJ), with the critical exponent 
P = — 1. Interestingly enough, the critical exponent /3 = 0, corresponding to a simple pole of the associated 
generating function, is otherwise known to arise in several simplified models, such as column-convex, convex, 
and directed polyominoes [H HI [5D] . 

As the foregoing discussion suggests, there is considerable interest in solving, exactly, probabilistically, or 
asymptotically, restricted models of self-avoiding walks and polygons. Beyond serving to develop informed 
conjectures regarding more complex models, this is relevant to areas such as statistical physics and the 
statistical mechanics of polymers |36j . For combinatorialists, we may observe that consideration of such 
models has served as a powerful incentive to develop new counting methods based on generating functions [U 
[3 [201 EH] i including transfer matrix methods and what is known as the "kernel method" . 

The present article focuses on a special type of self-avoiding polygons, the 3-sided prudent polygons (to 
be defined shortly - see Definition [I] in Section [2]) , when these are enumerated according to area. Roughly, 
a walk is prudent if it never takes a step towards an already occupied site and it is 2-, 3-, 4-sided if it has, 
respectively, 2, 3, or 4 allowed directions of growth; a prudent polygon is a prudent walk that is almost 
closed. For area n, we will obtain a precise asymptotic formula (Theorem [2] below) , 

(3) PA n ~ C(n) • 2" • n 9 , 

one that has several distinguishing features: (i) the critical exponent is the transcendental number g = log 2 3, 
in sharp contrast with previously known examples where it is invariably a "small" rational number; (ii) the 
multiplier C is no longer a constant, but a bounded quantity that oscillates around the value 0.10838 . . . and 
does so with a minute amplitude of 10~ 9 . The oscillations cannot be revealed by any standard numerical 
analysis of the counting sequence PA n , but such a phenomenon may well be present in other models, and, 
if so, it could change our whole view of the asymptotic behaviour of such models. 

Plan and results of the paper. Prudent walks are defined in Section [2j where we also introduce the 
2-sided, 3-sided, and 4-sided varieties. The enumeration of 2- and 3-sided walks and polygons by perimeter 
is the subject of insightful papers by Bousquet-Melou [5] and Schwerdtfeger [55] who obtained both exact 
generating function expressions and precise asymptotic results. In Subsections |2.2[^2.4| we provide the 
algebraic derivation of the corresponding area results: the enumeration of 2- and 3- sided polygons according 
to area is treated there; see Theorem[T]for our first main result. For completeness, we also derive a functional 
equation for the generating function of 4-sided prudent polygons (according to area), which parallels an 
incremental construction of [3J §6.5] - this functional equation suffices to determine the counting sequence in 
polynomial time. Section [3] dedicated to the asymptote analysis of the number of 3-sided prudent polygons, 
constitutes what we feel to be the main contribution of the paper. We start from a g-hypergeometric 
representation of the generating function of interest, PA(z), and proceed to analyse its singular structure: 
it is found that PA{z) has poles at a sequence of points that accumulate geometrically fast to |; then, 
the Mellin transform technology [T7] provides access to the asymptotic behaviour of PA(z), as z — > | in 
extended regions of the complex plane. Singularity analysis Ch.VI-VII] finally enables us to determine 
the asymptotic form of the coefficients PA n (see Theorem[2]) and even derive a complete asymptotic expansion 
(Theorem |3j. As already mentioned, the non-standard character of the asymptotic phenomena found is a 
distinctive feature. Section [4] concludes the paper with brief remarks relating to asymptotic methodology. 
In particular, experiments suggest that similar asymptotic phenomena are likely to be encountered in the 
enumeration of 4-sided prudent polygons. 

A preliminary announcement of the results of the present paper is the object of the communication pQ. 

2. Prudent walks and polygons 

One interesting sub-class of self-avoiding walks (SAWs) for which a number of exact solutions have been 
recently found are prudent walks. Introduced by Prea |30j . these are SAWs which never take a step towards an 
already occupied node. Exact solutions of prudent walks on a 2-dimensional square lattice were later studied 
by Duchi [llj and Bousquet-Melou [3], who were able to enumerate certain sub-classes. The enumeration of 
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FIGURE 1. Examples of (A) a two-sided prudent SAW; (B) a three-sided prudent SAW; 
(C) an (unrestricted) prudent SAW; (D) a prudent SAW leading to a prudent SAP. 



the corresponding class of polygons is due to Schwerdtfeger [32] • In this section, we first recall the definition 
of prudent walks and polygons, then summarize the known results of [3, 32 relative to their enumeration 
by length or perimeter; see Subsection |2.1[ where 2-, 3-, and 4-sided prudent walks are introduced. We 
then examine the enumeration of polygons according to area, in each of the three non-trivial cases. The 
case of 2-sided polygons is easy enough (Subsection 2.2). The main result of this section is Theorem [I] of 



Subsection |2.3[ which provides an explicit generating function for 3-sided polygons - it is on this expression 
that our subsequent asymptotic treatment is entirely based. In the case of 4-sided (i.e., "general") polygons, 
we derive in Subsection |2.4| a system of functional equations that determines the generating function and 
amounts to a polynomial-time algorithm for the generation of the counting sequence. 

2.1. Main definitions and results. We use the same classification scheme as the authors of [3j [32] - By 
definition, the endpoint of every prudent walk always lies on the boundary of the smallest lattice rectangle 
which contains the entire walk, referred to here as the bounding box or just box. This property leads to a 
natural classification of prudent walks (see Figure [I]) . 

Definition 1. Let u be a prudent walk of length n, and let u>i be the prudent walk comprising the first i 
steps of L). Let bi be the bounding box of u>i. Then uj is 1-sided if uji ends on the north side of bi for each 
i = 0, 1, . . . , n; 2-sided if each u>i ends on the north or east sides of bi; 3-sided if each u>i ends on the north, 
east or west sides ofbi (with one caveat, described below); and 4-sided (or unrestricted,) if each uji may end 
on any side of bi . 

Remark. The issue with 3-sided prudent walks is encapsulated by the walk (0,0) —> (1,0) — > (1,-1) — > 
(0, —1). If one draws the walk's box after each (discrete) step, then it is clear that the walk always ends on 
the north, east or west sides, seemingly fulfilling the 3-sided requirement. However, if the walk is taken to 
be continuous, then along the step (1, —1) — > (0, —1), the endpoint is only on the south side. In general this 
occurs when a walk steps from the south-east corner of its box to the south-west corner (or vice versa) when 
the box has width one and non-zero height. Allowing such walks forces us to account for structures like 
those in Figure [2] while this is certainly possible, it complicates a number of rational terms and contributes 
little to the asymptotic behaviour of the model. For this reason we follow the examples of Bousquet-Melou 
and Schwerdtfeger [3j [32] and exclude these cases. 

An equivalent definition of the 3-sided walks considered here is to forbid two types of steps: when the 
box has non-zero width, a south step may not be followed by a west (resp. east) step when the walk is on 
the east (resp. west) side of the box. 



Walks (length). The enumeration of 1-sided prudent walks (also known as partially directed walks) is 
straightforward, and the generating function for such walks is rational - we will not discuss these any 
further. Duchi [TT] successfully found the generating function for 2-sided walks, showing it to be algebraic. 
Bousquet-Melou [3] solved the problem of 3-sided prudent walks, finding the generating function to be non- 
D-finite. For the unrestricted case, functional equations were found by both Duchi and Bousquet-Melou, but 
at present these equations remain unsolved. 
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Figure 2. Examples of the (A) prudent walks and (B) prudent polygons which we exclude 
from the definition of 3-sided. 



The dominant singularity of the generating functions for 2-sided and 3-sided prudent walks is a simple 
pole at p = 0.4030317168..., the smallest real root of 1 — 2x - 2x 2 + 2x 3 . Dethridge and Guttmann 
conjecture that the same is true for unrestricted prudent walks, based on a computer-generated series of 100 
terms. They also conjectured that the generating function for unrestricted prudent walks is non-holonomic[j] 
(non-D-finite). Here is a summary of known results, with the estimates tagged with a question mark (?) 
being conjectural ones: 



(4) 



Walks 


Generating function 




Asymptotic number 


References 


2-sided 


algebraic 


t2 


p~ n , p~ l ~ 2.481 


Duchi QT], Bousquet-Melou [3] 


3-sided 


non-holonomic 




p~ n , p- 1 ~ 2.481 


Bousquet-Melou [3] 


4-sided 


functional equation 


K4 


p~ n , p" 1 - 2.481 (?) 


Dethridge & Guttmann [jjj. 



The values of the multipliers, after [9j, are K2 = 2.51..., K3 = 6.33... and (estimated) K4 s» 16.12. 

Prudent polygons (perimeter). Self-avoiding polygons (SAPs) are self-avoiding walks which end at a 
node adjacent to their starting point (excluding walks of a single step). If the walk has length n — 1 then 
the polygon is said to have perimeter n. 

Definition 2. A prudent self-avoiding polygon (prudent polygon) is a SAP for which the underlying SAW 
is prudent. In the same way, a prudent polygon is 1-sided (resp. 2-sided, etc.) if its underlying prudent walk 
is 1-sided (resp. 2-sided, etc.). 

A 1-sided prudent polygon starting at (0, 0) must end at (0, 1), thus consisting only of a single row of cells 
and having a rational generating function. The enumeration of 2- and 3-sided prudent polygons by perimeter 
has been addressed by Schwerdtfeger |32j . The non-trivial 2-sided prudent polygons are essentially inverted 
bargraphs [3T] , and so the 2-sided case has an algebraic generating function. Schwerdtfeger finds the 3-sided 
prudent polygons to have a non-D-finite generating function. 

If PPW(z) = J2 n >o 

pit* z n is defined to be the half-perimeter generating function for fc-sided prudent 

polygons (so p$ is the number of fc-sided prudent polygons with perimeter 2n), then the following holds [32] : 
(i) the dominant singularity of PP^(z) is a square root singularity at a = 0.2955977..., the unique real 
root of 1 — 3a; — x 2 — x 3 . So pi? ~ A2(7~ ,l n~ 3 / 2 as n —> 00, where A2 is a constant; (ii) the dominant 
sm gularity of PP^ (z) is a square root singularity at r = 0.24413127..., where r is the unique real root of 
t 5 + 6t 3 — 4t 2 + 17r — 4. So p^ ~ A3T~"?i -3 / 2 as n — > 00, where A3 is a constant. Schwerdtfeger has 
furthermore classified 4-sided prudent polygons in such a way as to allow for functional equations in the 
generating functions to be written. Unfortunately no one has thus far been able to obtain a solution from 
said equations. We again present a summary of known results. 



A function is said to be holonomic or D-finite if it is the solution to a linear differential equation with polynomial 
coefficients; see [33] Ch. 6] and [20] Ap. B.4]. 
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(5) 



Polygons (perimeter) 


Generating function 




Asymptotic 


number 


References 


2-sided 


algebraic 


A 2 


-n -3/2 
g n , a 


_1 ~ 3.382 


Schwerdtfeger [32] 


3-sided 


non-holoriomic 


A 3 


-n -3/2 
r n , t 


_1 ~ 4.096 


Schwerdtfeger [32] 


4-sided 


functional equation 


A 4 


— n —5 

v n , V 


_1 « 4.415 (?) 


Dethridge et al. [S] 



The empirical estimates regarding 4-prudent polygons are taken from [5]. They are somewhat imprecise, 
and it is suspected that the critical exponent satisfies 5 — —3.5 ± 0.1, with S = —7/2 a compatible value. 

Prudent polygons (area). The focus of this paper is on the enumeration of prudent polygons by area, 
rather than perimeter. The constructions we use here are essentially the same as Schwerdtfeger's [35]; the 
resulting functional equations and their solutions, however, turn out to be quite different, as will be revealed 
by the peculiar singularity structure of the generating functions and the non-trival asymptotic form of the 
coefficients. We have modified Schwerdtfeger's construction for 4-sided prudent polygons slightly to allow 



for an easier conversion into a recursive form (see Subsection 2.4 1. 

We will denote the area generating function for fc-sided prudent polygons by PA^ k \q) = XL>1 PA n q n . 
For 3- and 4-sided prudent polygons, it is necessary to measure more than just the area - in these cases, 
additional catalytic variables will be used (see [3] for a more thorough explanation). 

2.2. Enumeration of 2-sided polygons by area. The non-trivial 2-sided prudent polygons can be con- 
structed from bargraphs. Let B(q) = J2 n >i t>nq n be the area generating function for these objects. The area 
generating function for bargraphs, B(q), is 

and so b n = 2 n ~ 1 for n > 1. (Bargraphs are a graphical representation of integer compositions.) 
Proposition 1. The area generating function for 2-sided prudent polygons is 

PA^(q) = -\ + - 2 ^, 
1 — 2q 1 — q 

(21 

and so the number of such polygons is PAn = 2™ + 2 for n > 1. 

Proof. A 2-sided prudent polygon must end at either (0, 1) or (1,0). Reflection in the line y = x will not 
invalidate the 2-sided property, so it is sufficient to enumerate those polygons ending at (1,0) and then 
multiply the result by two. 

The underlying 2-sided prudent walk cannot step above the line y = 1, nor to any point (x, y) where 
x, y < 0. So any polygon beginning with a west step must be a single row of cells to the left of the y-axis. 
The generating functions for these polygons is then -^L- , 

A polygon starting with a south or east step must remain on the east side of its box until it reaches the 
line y = 1, at which point it has no choice but to take west steps back to the y-axis. It can hence be viewed 
as an upside-down bargraph with north-west corner (0, 1). The area generating function for these objects is 

Adding these two possibilities together and doubling gives the result. □ 

2.3. Enumeration of 3-sided polygons by area. When constructing 3-sided prudent polygons, we will 
use a single catalytic variable which measures width. To do so we will need to measure bargraphs by width. 
Let B(q 1 u) = X) n >i Si>i ^n,i9 n, w l be the area- width generating function for bargraphs (so b n ^ is the number 
of bargraphs with area n and width i). 

The area-width generating function for bargraphs, B(q,u), satisfies the equation 

(6) B (q,u) = ^- + ^-B(q,u), 

which is obtained by successively adding columns. Accordingly, by solving the functional equation, we obtain 

7-./ \ 1 U 



L 



+ 



FIGURE 3. The decomposition used to construct 3-sided prudent polygons. 



and so b n ^ = for n,i> 1. (Clearly, 6 ni counts compositions of n into i summands.) 

Let W(q,u) — J2 n >i Si>i w n,iq n u l be the area-width generating function for 3-sided prudent polygons 
which end at (—1, 0) in a counter-clockwise direction. As we will see, this is the most complex type of 3-sided 
prudent polygon; everything else is either a reflection of this or can be constructed from something simpler. 

Lemma 1. The area-width generating junction for 3- sided prudent polygons ending at (—1,0) in a counter- 
clockwise direction, W(q,u), satisfies the functional equation 

(7) W(q, u) = qu(l + B(q, «)) + -*—(W(q, u) - W(q, qu)) + qu{l + B(q, u))W(q, qu) 

1-q 

Proof. The underlying prudent walk cannot step to any point (x, y) with x,y < 0, nor to any point with 
x < — 1. It must approach the final node (—1,0) from above. So the only time the endpoint can be on the 
west side of the box and not the north or south is when the walk is stepping south along the line x = — 1. 
So prior to reaching the line x = — 1, the walk must in fact be 2-sided. Note that the north-west corner of 
the box must be a part of of the polygon. 

If the walk stays on or below the line y = 1, then (as will be seen in Proposition [I]) , it either reaches the 
point (0, 1) with a single north step, or by forming an upside-down bargraph. This must then be followed by 
a west step to (— 1, 1), then a south step. This will form either a single square or a bargraph with a single 
square attached to the north-west corner, giving the first term on the right-hand side of d7|). 

Since the north-west corner of the box of any of these polygons is part of the polygon, it is valid to add 
a row of cells to the top of an existing polygon (so that the west sides line up). This can be done to any 
polygon. If the new row is not longer than the width of the existing polygon we obtain the term 

n>l i>l k=l n>l i>l " 

giving the second term in the right-hand side of (J7|. 

Note. For the remainder of this subsection, we will omit unwieldy double or triple sums like the one above, 
and instead give recursive relations only in terms of the generating functions. 

Instead, the new row may be longer than the width of the existing polygon. In this case, as the walk steps 
east along this new row, it will reach a point at which there are no occupied nodes south of its position, and 
it will hence be able to step south in a prudent fashion. It must then remain on the east side of the box 
until reaching the north side, at which point it steps west to x = —1 and then south to the endpoint. This 
effectively means we have added a row of length equal to the width +1, and then (possibly) an arbitrary 
bargraph. So we obtain 

quW(q,qu)(l + B(q,u)) 

which gives the final term in the right-hand side of Q. □ 

Lemma 2. The area-width generating function for 3-sided prudent polygons ending at (—1,0) in a counter- 
clockwise direction is 

oo m— 1 

W(q,u) = Y F{q,q m u) J] G(q,q k u), 
m=0 k=Q 
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where 

qu(l-q) 2 -q(l-q-u + qu-q 2 u) 

= (1-2 9 )(1 -q-quY G ^ = { l-2q)(l-q-qu) " 

Proof. Substituting BW(q, u) — 1 _^" gM into |7| and rearranging gives 

(8) W(q,u) = F(q,u) + G(q,u)W(q,qu) 
Substituting u i-» uq gives 

(9) W(q, qu) = F(q, qu) + G(q, qu)W{q, q 2 u) 
and combining these yields 

(10) W(q, u) = F(q, u) + F(q, qu)G(q, u) + G(q, u)G{q, qu)W(q, q 2 u). 
Repeating for u i— > q 2 u, q 3 u, q M u will give 

M m—1 M 

(11) W(q, U)=J2 F (l> <T«) II G ^> ?*«) + II G ^ <? n u)W{q, q M+1 u). 

m=0 fc=0 m=0 

We now seek to take M — > oo. To obtain the result stated in the Lemma, it is necessary to show that 

M ro-1 
^F(q,g™ W ) [] G(q,q k u) 
m=0 fc=0 

converges, and 

M 

Y[ G{q,q m u)W{q,q M+1 u) ->■ 

m=0 

as M — > oo (both considered as power series in q and u). 

Both F and G are bivariate power series in q and u. We have that 

F(q, u) = qu + q 2 {u + u 2 ) + q 3 {2u + 2u 2 + u 3 ) + 0{q 4 ) 

G{q, u) = q(-l + u) + q 2 {-2 + u + u 2 ) + g 3 (-4 + 2u + 2u 2 + u 3 ) + 0{q A ) 
It follows that F(q, q m u) = 0(q m+1 ) and G(q, q k u) = O(q) for all m, k > 0. So we have that 

m— 1 

F(q,q m u) J] G(q,q k u) = 0(q 2m+1 ) 

k=0 



So considered as a power series in q and u, the first term in the right-hand side of (111 does converge to a 
fixed power series as M —¥ oo. 

By the same argument, we have that 

M 

[] G(q,q m u)^0 
m— 

as M — > oo. So it suffices to show that W(q, q M+1 u) converges to a fixed power series. But now every term 
in the series W(q, u) has at least one factor of u (since every polygon has positive width), so it immediately 
follows that W(q,q M+1 u) ->• as M -> oo. 



So both terms in ( 11 ) behave as required as M — ¥ oo, and the result follows. □ 



Theorem 1. The area generating function for 3-sided prudent polygons is 

9/7 3 n n) 2 J2. < l\m 2m m ~} 1 „ n k i rl k+l n k+2 

pA(3)( \ _ ~ z 1 I 1 ~ 1) \ L ) g TT 1 9 q +q ~ q 

W (1-2(7)2 ^ (1 - 2g) w (l - q - q m+1 ) ^ 1 - 4 - g fc+1 

2g(3 - lOg + 9g 2 - q 3 ) 
+ (l-2q) 2 (l-q) 
= 6q + 10q 2 + 20q 3 + 42q 4 + 92q 5 + 204g 6 + 454q 7 + lOlOg 8 + 2242 9 9 + 4962g 10 ... 



Proof. A 3-sided prudent polygon must end at (—1,0), (0,1) or (1,0), in either a clockwise or counter- 
clockwise direction. Setting u — 1 in W(q,u) gives the area generating function 

no) w( u - q3(1 ~ q)2 V (-1) V m "fi 1 1 - g q k + g fc+1 - q k+2 | g(i - q? 

[ ' W ' j (l-2 9 ) 2 ^ (l-2g)" l (l-g-?"» +1 ) JUL 1 - q - q k + l (l-2g) 2 ' 

A clockwise polygon ending at (—1,0) can only be a single column, which has generating function 
(13) 



1-3 

A counter-clockwise polygon ending at (0, 1) cannot step left of the y-axis or above the line y = 1. While 
it is below this line, it must remain on the east side of its box, and upon reaching the line y = 1, it must 
step west to the y-axis. It must therefore be a bargraph, with generating function 

<») rhi- 

A reflection in the y-axis converts a polygon ending at (—1,0) to one ending at (1,0) in the opposite 
direction, and reverses the direction of a polygon ending at (0, 1). So adding together and doubling (12 1, 
(13) and (14) will cover all possibilities, and gives the stated result. □ 

2.4. Enumeration of 4-sided polygons by area. This case is included for completeness, as the results 
are not needed in our subsequent asymptotic analysis. A 4-sided prudent polygon may end at any of 
(0, 1), (1, 0), (0, — 1), (— 1, 0) in either a clockwise or counter-clockwise direction. Reflection and rotation 
leads to an 8-fold symmetry, so it suffices to count only those ending at (—1,0) in a counter-clockwise 
direction. We modify Schwerdtfeger's sub-classification slightly. 

Let X{q,u,v) = J2 n >i Si>i J2j>i x n,i,jq n u l vi be the generating function for those polygons (class X) 
for which removing the top row does not change the width or leave two or more disconnected pieces, with q 
measuring area, u measuring width and v measuring height. 

Let Y(q, u, v) — J2 n >i Si>i ^2j>i Vn.i^q 11 " 1 ^^ be the generating function for the unit square plus those 
polygons (class y) not in X for which removing the rightmost column does not change the height or leave 
two or more disconnected pieces, with q measuring area, u measuring height and v measuring width. 

Let Z{q,u,v) = J2 n >i Si>i Sj>i z n.i.j < f' ulv ^ be the generating function for the polygons (class Z) not 
in X or y, with q measuring area, u measuring width —1 and v measuring height. 



Proposition 2. The generating functions X(q, u, v), Y(q, u, v) and Z(q, u, v) satisfy the functional equations 
(15) 



qv qv 
X(q, u, v) = j— [X(q, u, v) - X(q, qu, v)] + — - [Y(q, v, u) - Y(q, v, qu)} 



q 

quv 



i _ [Z(q,u,v) - qZ(q,qu,v)] 
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qv qv 
Y(q,u,v) = quv + [Y(q,u,v) - Y(q,qu,v)} + [Z(q,v,u) - Z(q,v,qu)} 

(16) l-q 1-q 

+quv[X(q, qv, u) + Y(q, u, qv) + qvZ(q, qv, u)] 

qv 

(17) Z(q, u, v) = [Z(q, u, v) - Z(q, qu, v)\ + qvY(q, qv, u) + quvZ(q, u, qv) 

The generating function for 4-sided prudent polygons is then given by 

PA^(q) = 8[X(q, 1, 1) + Y(q, 1, 1) + Z(q, 1, 1)] 

= 8q + 24q 2 + 80q 3 + 248g 4 + 736q 6 + 2120g 7 + 5960g 8 + 16464g 9 + 44808q 10 + . . . 

Proof. As with the 3-sided polygons in Lemma [lj the walk cannot visit any point [x, y) with x, y < or 
with x < — 1. The walk must approach (—1,0) from above, and must do so immediately upon reaching the 
line x = — 1. So every polygon contains the north-west corner of its box. As in the 3-sided case, this leads 
to a construction involving adding rows to the top of existing polygons. 



X 



X+Y+Z 





Figure 4. The decompositions used to construct 4-sided prudent polygons in X, y, Z 
(from top to bottom). 

By definition, a polygon in X of width i can be constructed by adding a row of length < i to the top of 
any polygon of width i. Adding a row to a polygon in X gives 

EEE zn.^w • v q k = qv 12 • 

n>l i>l j>l k=l n>l i>l j>1 ^ 



Z gives the rest of (151 



which is the first term in the right-hand side of (15). Performing similar operations for polygons in y and 



Note. Again, for the remainder of this subsection we give recursive relations purely in terms of the generating 
functions. 

Polygons not in X must also contain the north-east corner of their box. This leads to another construction 
involving adding columns to the right-hand side of existing polygons. To obtain a polygon in y of height i, 
a new column of height < i should be added to a polygon of height i which contains the north-east corner 
of its box. So adding a column to a y polygon gives 

1 V vr \ 1 y vi \ 
-Y(q,u,v) - - Y(q,qu,v) 



1-0 ' ' 1 



which is the second term in the right-hand side of (16 1. Performing a similar operation for Z polygons gives 



the third term in (16 1. 

Adding a new column to a polygon in X containing its north-east corner can be viewed as adding a 
sequence of rows on top of one another, and so if the new column has height > 2 then the resulting polygon 
is actually in X. If the new column has height one, however, the resulting polygon is in y. Isolating those 
polygons in X which contain their north-east corner is difficult; however, we can perform an equivalent 
construction by adding a row of length i + 1 to any polygon of width i. Doing so to a polygon in X gives 

quvX(q, qv, u) 



and combining this with the same for y and Z gives the fourth term in (16). The quv term is the unit 
square. 

Polygons in Z also contain the south-east corner of their box. In a similar fashion to the constructions 
for y and Z, we can add a new row to the bottom of a polygon containing its south-east corner. To do so 

9 



to a polygon in Z of width i + 1 (remember u measures width —1) requires a new row of width < i, so we 
obtain 

Z(q,u,v) - - Z{q,qu,v) 



1 



which is the first term in the right-hand side of (17 1 



Adding a new row to the bottom of something in X (containing its south-east corner) will give back 
something in X ', which will have been constructed by an alternate method described above. Adding a new 
row of length > 2 to a polygon in y will result in another polygon in y, which will also be constructible via 
alternate means. So we are left only with the possibility of adding a row of length one to the bottom of a 
polygon in y. This is analogous to the above description of adding a column of height one to the right of a 
polygon in X; we now proceed by adding a column of height i + 1 to a polygon in y or Z of height i. Doing 



so gives the final two terms in (17). □ 



3. Asymptotics 

For most lattice object problems, finding and solving the functional equation(s) is the difficult part. Once 
a generating function has been found, the dominant singularity is often quite obvious, and so the asymptotic 
form of the coefficients can be easily described. The problem of 3-sided prudent polygons, however, turns out 
to be rather the opposite. The functional equation ([T]) was not terribly difficult to obtain, and its solution 
is relatively simple - it only comprises a sum of products of rational functions of q. 

The asymptotic behaviour of this model, on the other hand, is considerably more complex than any model 
we have seen before. The dominant singularity at q = 1/2 is not even apparent from the representation of 
Theorem [I] As we shall see, there is in fact an accumulation of poles of the generating function^] PA (z) 
towards q = 1/2. Accordingly, the nature of the dominant singularity at q — 1/2 is rather unusual: a singular 
expansion as q approaches 1/2 can be determined, but it involves periodic fluctuations, a strong divergence 
from the standard simple type Z a (\og Z)@ , where Z := 1 — z/p, with p (here equal to 1/2) the dominant 
singularity of the generating function under consideration. This is revealed by a Mellin analysis of PA(z) 
near its singularity, and the periodic fluctuations, which appear to be in a logarithmic scale, eventually echo 
the geometric speed with which poles accumulate at 1/2. Then, thanks to a suitable extension to the complex 
plane, the singular expansion can be transfered to coefficients by the method known as singularity analysis [301 
Ch. VI]. The net result is, for the coefficients PA n , an asymptotic form that involves a standard element 2 n n 9 , 
but multiplied by a periodic function in log 2 n. The presence of oscillations, the transcendental character of 
the exponent g — log 2 3, and the minute amplitude of these oscillations, about 10~ 9 , are noteworthy features 
of this asymptotic problem. 

Theorem 2. The number PA n = PA„ of 3-sided prudent polygons of area n satisfies the estimate 
(18) PA n = [k + K{log 2 n)] 2 n -n 9 + 0(2 n -n^Hogn) , n -> oo, 

where the critical exponent is 



and the "principal" constant is 
(19) K ° = 91oe2sin(™)r(o + l) H 



g = log 2 3 = 1.58496 

(l-§2-')(l-§2-*) 



9 log 2 sin(7r 5 )r( g + 1) 11 (1 - |2-i)a 



0.10838 42946. 



The function n(u) is a smooth periodic function of u, with period 1, mean value zero, and amplitude = 
1.54623 • 10 -9 , which is determined by its Fourier series representation: 

, s ?,k7T„ •„ sin(7ro) r(l + g) 

k(u) = \ K fe e il , with n k = k 



fc J^ ni — - sin(7r 9 + 2z^ 2 /log2) r(l + 3 + 2ifc7r/log2)' 



2 Throughout this section only dedicated to 3-sided prudent polygons, we omit redundant superscripts and let PA n and 



PA(q) represent, respectively, what was denoted by PA^ and PA^(q) in Section^ 
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Quantity 



atq = 1/2 



reference 



l-9 + i 



1 - i 
q 2 



7 = i 
C(q) 



l-q + q 2 
log v 



J 1/9 

2g(3 - lOg + 9g 2 - g 3 ) 

(l-g)(l-2g)2 
2g(l- g r 

(l-2g) 2 



3 
2 
1 
3 

log 2 (3/2) 
1 

" 4(1 - 2g) 2 
1 

" 4(1 - 2g) 2 



Eq. (231 



Eq. (231 



Eq. (411 



Eq. (501 



Eq. (211 and (26 1 



Eq. (211 and (25 1 



FIGURE 5 . A table of some of the recurring quantities of Section [sj their reduction at 
and the relevant equations in the text. 



1/2 



The proof of the theorem occupies the next subsections, whose organization reflects the informal descrip- 
tion given above. We shall then discuss the fine structure of subdominant terms in the asymptotic expansion 
of PA n ; cf Theorem [3] Some quantities that appear repeatedly throughout this section are tabulated in 
Figure [5] for convenience. 

3.1. Resummations. We start with a minor reorganization of the formula provided by Theorem [TJ com- 
pletion of the finite products that appear there leads to the equivalent q-hypergeometric form 



(20) PA(q)=C(q)+A(q)-Q(l;q)-J2^ 



(-l)V 



^(l-2g)» Q(q n ;q)- 
Here and throughout this section, the notations are 

(21) C(a) - M3-Wq + % 2 -<l 3 ) A( ) . gg(l - g? 

and 

(vz;q) c 



(22) Q(z;q) :=Q(z;q;u(q),v(q)) , where Q(z; q; u, v) 



(quz^q)^ 
with 

(23) U(q)=« V(q)= 1 -^±f. 

1-q 1-q 

In the definition of Q, the notation (a;; q) n represents the usual g-Pochhammer symbol: 

(x; q) n = (1 - x){\ -qx)---{l- xq"' 1 ). 

Lemma 3. The function PA(q) is analytic in the open disc \q\ < s/2 — 1, where it admits the convergent 
g-hypergeometric representation 

(24) PA(q) = C(q) + A(q)p^ f>l)« ■ ^^C^k , 



with A(q),C(q),u = u(q),v = v(q) rational functions given by (|21[) and (23) 



Proof. The (easy) proof reduces to determining sufficient analyticity regions for the various components of 



the basic formula (20), some of the expansions being also of later use. First, the functions A(q) and C(q) 



are meromorphic for |g| < 1, with only a pole at q = 1/2. They can be expanded about the point q = 1/2 



to give 

1 1 1 1 - 2q 

( 5) ~ 4(1 - 2q) 2 + 4(1 - 2q) ~ 4 4~ 

The function Q(l; q) is analytic for \q\ < 1 except at the points for which (uq; q)^ = 0, that is, the points 
a for which 1 — a — <r n = for n > 2. The smallest of these (in modulus) is ip = (V5- l)/2 = 0.618034..., 
a root of 1 — q — q 2 . So Q(l;g) is certainly analytic at q — 1/2; the constant term in its expansion about 
q = 1/2 is 

« 1;1/2 > = §f^ = -»' 18109m -- 

In similar fashion, 1/Q(z; q) is bivariate analytic at points {z, g) for which \q\ < 1, except when (vz; q , ) 00 = 
0. This occurs at points (zj,q) where Zj := for j > 0. In particular, for \q\ < 9, wher^~~ 

(27) 6» = 0.56984 := the unique real root of 1 - 2x + x 2 - x 3 , 

we have \z<j\ > 0, hence \zj\ > 9, for all j > 0. So, X/Q{z;q) is analytic in the region {(z,q) : \z\, \q\ < 9}. 
Thus, for all n > 1, the functions 1/Q(q n ; q) are all analytic and uniformly bounded by a fixed constant, for 
\q\ < r , where r is any positive number such that r < 9. 

From these considerations, it follows that the central infinite sum that figures in (20) is, when \q\ < r%, 
dominated in modulus by a positive multiple of the series 

(28) E(i_2 ri )«' 

provided that r\ < 9 and r\/{l — 2ri) < 1. Any positive r\ satisfying r\ < \/2 — 1 is then admissible. 
In that case, for \q\ < r\, the central sum is a normally convergent sum of analytic functions; hence, it is 
analytic. □ 

The radius of analyticity of PA(q) is in fact 1/2. In order to obtain larger regions of analyticity, one needs 
to improve on the reasoning underlying the derivation of ( 28 1 . This will result from a transformation of the 
central infinite sum in ( 20 ) , namely, 

,2n 1 



(29) %):=£(-!)"- 



(l-2<?)« Q(q n ;q)' 

Only the bound 1/Q(q n ; q) = 0(1) was used in the proof of Lemma[3j but we have, for instance, 1/Q(q n ; q) = 
1 + 0(q n ), as n — > oo, and a complete expansion exists. Indeed, since 1/Q is bivariate analytic in |z|, \q\ < 9, 
its z-expansion at the origin is of the form 

(30) _^ = 1 + 5>(ffK . 



In particular, at z = q n , we have 



(3D 



q vn . 



U>1 



Now, consider the effect of an individual term d v (q) (instead of \/Q(q n ;q)) on the sum (29 1. The identity 

n 2n rt vJr2 

^ E(-i) w (r^^ = - 1 - 2 ^ +r , a 



3 The function v(q) = 1 + g 2 /(l — q), having nonnegative Taylor coefficients, satisfies \v{q)\ < v(\q\), for \q\ < 1; thus, 
| | > Also, l/v(x) decreases from 1 to for x £ [0, 1]. Hence, with 9 the real root of l/v(8) = 8, it follows that 

\zo\ > 8 as soon as |g| < 8. 
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provides an analytic form for the sum on the left, as long as q is not a pole of the right-hand side. Proceeding 



formally, we then get, with (31 1 and (32 ), upon exchanging summations in the definition (29 ) of S(q), a form 
of PA(q) that involves infinitely many meromorphic elements of the form 1/(1 — 2q + q u+2 ). 

We shall detail validity conditions for the resulting expansion; see (34) below. What matters 

2q + q»+ 2 ) 



from ( 32 ) , is the location of poles of the rational functions (1 
(33) 



as seen 

for v > 1. Define the quantities 



Cfe := the root in [0, 1] of 1 - 2a; 



r /c+2 



0. 



We have 



Co = i; Ci 



\/5-l 



0.618, C2 = 0.543, Ca = 0.518, 



and Cfe — > \ as k increases. The location of the complex roots of 1 — 2x + x k+2 = is discussed at length 
in [5D1 Ex. V.4, p 308], as it is related to the analysis of longest runs in binary strings: a consequence of 
the principle of the argument (or Rouche's Theorem) is that, apart from the positive real root Cfc) all other 
complex roots lie outside the disc \z\ < |. The statement below builds upon this discussion and provides 
an extended analyticity region for PA(q) as well as a justification of the validity of the expansion resulting 
from (31) and (32), which is crucial to subsequent developments. 



Lemma 4. The generating function PA(q) is analytic at all points of the slit disc 

Vo:={q : \q\ < [± } ■ 

For geDo, the function PA(q) admits the analytic representation 



(34) 



where 



PA(q) = C(q) - A(q) 



(qu;q) c 



q 



(l-q)" 



d v {q) 



l-2q + q 



v+2 



du(q) 



I 1 = r z v, (quz;q)c 

Q{z;q) (vz;q) x 



In the disc \z\ < 



55 

100 



punctured at \, the function PA(q) is meromorphic with simple poles at the points 



CsjCsj'-'j with Cfe as defined in (33). Consequently, the function PA(q) is non-holonomic, and, in particular, 
transcendental. 

Proof. The starting point, noted in the proof of lemma [3j is that fact that 1/Q(z;q) is bivariate analytic at 
all points (z,q) such that \z\, \q\ < 8, where 8 = 0.56984 is specified in (27). Cauchy's coefficient formula, 



d v {q) 



dz 



is applicable for any 81 such that < 8 X < 8. Let us set 81 = j^. Then, since 1/Q(z;q) is analytic, hence 
continuous, hence bounded, for \z\ < 8\ and \q\ < 8\, trivial bounds applied to the Cauchy integral yield 



(35) 

for some absolute constant C > 0. 



K(q)\<c- 



Consider the double sum resulting from the substitution of ( |31[ ) into ( 29 ) 



n>l 



(1 - 2q)' 



If we constrain q to be small, say \q\ < j^, we see from (35) that the double sum is absolutely convergent. 
Hence, the form (34) is justified for such small values of q. We can then proceed by analytic continuation 



from the right-hand side of ( 34 ) . The bound ( 35 ) grants us the fact that the sum that appears there is indeed 



analytic in T>q. The statements, relative to the analyticity domain and the alternative expansion (34) follow. 
Finally, since the value \ corresponds to an accumulation of poles, the function PA(q) is non-holonomic 
(see, e.g., [16] for context). □ 
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As an immediate consequence of the dominant singularity being at ^, the coefficients PA n must obey a 
weak asymptotic law of the form 

PA n ^2 n 6(n) 1 where limsup0(n) 1/ri = 1, 

n— >oo 

that is, 0(n) is a (currently unknown) subexponential factor. 

More precise information requires a better characterization of the behaviour of S(q), as q approaches the 
dominant singularity | . This itself requires a better understanding of the coefficients d v (q) . To this end, we 
state a general and easy lemma about the coefficients of quotients of g- factorials. 

Lemma 5. Let a be a fixed complex number satisfying \a\ < 1 and let q satisfy \q — -|| < j^. 
for v > 1 



One has, 



(36) 



(az;q) c 



E 



(aq ] ;q)c 



(z;q)oo (q;q)oo f^ (q J ;q) 3 

Proof. The function h(z) :— (az; q) oc /(z; q),^ has simple poles at the points Zj := q~ J , for j > 0. We have 



h(z) 



ej(a;q) 



1-zqi' 11 (q j ;q)j{q;q)oo 

The usual expansion of coefficients of meromorphic functions [20l Th. IV. 10, p. 258] immediately implies a 
terminating form for any J G Z> : 

(37) L .. - / ,',•!"• 'j "/ - ""'./'• 

where we may adopt Rj 

The last estimate ( 37 1 corresponds to an evaluation by residues of the Cauchy integral representation of 
coefficients, 



[z"]h(z)=^2 ej {a;q)qi v + 0(R 

3=0 



[z v )h{z) 



1 

2i7T 



h(z) 



\x\=Rj 



dz 

7 v+l ' 



Now, let J tend to infinity. The quantity Rj lies approximately midway between two consecutive poles, 
q~ J and q~ J_1 : and it can be verified elementarily that, throughout \z\ = Rj, the function h(z) remains 
bounded in modulus by an absolute constant (this requires the condition \a\ < 1). It then follows that we 
can let J tend to infinity in (37 1. For v > 1, the coefficient integral taken along \z\ — Rj tends to 0, so that, 
in the limit, the exact representation (36) results. □ 



The formula (36) is equivalent to the partial fraction expansion (Mittag-Leffler expansion; see [22] §7.10]) 
of the function h(z), which is meromorphic in the whole complex plane: 



(38) 



(az;q)c 

(z;q)c* 



l + 



(q;q) 



(aq J ;g)oo zq> 

»^ (q~ j \q)j ~L-zqi 



E 



(The condition \a\ < 1 ensures the convergence of this expansion.) As observed by Christian Krattenthaler 
(private communication, June 2010), this last identity is itself alternatively deducible from the q-Gaufi 
identitj^] 



2^1 



A,B 



G 



C ;q ' AB 



upon noticing that 



h(z) 



(q;g)oo 

(1 - z)(q\q) c 



(C/A;q) 00 (C/B;q) c 
(C-^UiC/iAB^q), 

q/a,z 



■ 291 



qz 



;q,a 



For notations, see Gasper and Rahman's reference text 121 1 : page 3 (definition of <f>) and Eq. (1.5.1), page 10 (g-Gaufi 
summation). 
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A direct consequence of Lemma [5] is an expression for the coefficients d v (q) — [z v ]Q{z; q) 1 , with Q(z; q) 
denned by p2l: 



(39) 

To see this, set 



E 



(quv q 1 : 



q)c 



(vq J Y 



v>\. 



quv 



1-q + q 2 



and replace z by zv in the definition of h(z). Note that at q = 1/2, we have u = 1, v = 3/2, a — 1/3, so 
that, for qfa 1/2, we expect d v (q) to grow roughly like (3/2)". 
Summarizing the results obtained so far, we state: 



Proposition 3. The generating function of 3-sided prudent polygons satisfies the identity 

'{aq- j ;q) 1 v v q^ +1 > 



(40) 

where 
(41) 



PA(q) = D(q) - q z A{q) 



(a;q)oo(v;q)oo 

u=l j=Q 



(?; q) 00 {av;q) c 



{q- ] ;q)j 



i-q 



i-q 



D(q) = C(q) 



1 - 2q + q»+ 2 
(v,q)oo 



A{q)- 



(l-q) 2 {av;q) 



and A(q),C(q) are rational functions defined in Equation (21). 



Proof. The identity is a direct consequence of the formula (39) for d v (q) and of the expression for PA(q) 
in ( 34 ) , using the equivalence av = qu and the simple reorganization 

(ag _J ;g)oo = (aq~ J ;q)j ■ (a;g)oo- 

Previous developments imply that the identity (40) is, in particular, valid in the real interval (0, |). The 
trivial equality 

U2) {aq~ 3 \q)j = (a-q)(a-q 2 )---{a-q J ) 

1 > (g-i;^- (l- q )(l- q 2)...(l- q j) 

then shows that the expression on the right-hand side indeed represents a bona fide formal power series in q, 
since the q- valuation of the general term of the double sum in (40) increases with both j and v. □ 

The formula (40) of Proposition [3] will serve as the starting point of the asymptotic analysis of PA(q) as 
q — > 1/2 in the next subsection. Given the discussion of the analyticity of the various components in the proof 
of Lemma [3j the task essentially reduces to estimating the double sum in a suitable complex neighbourhood 
oiq= 1/2. 



3.2. Mellin analysis. Let T{q) be the double sum that appears in the expression (40) of PA(q). We shall 
take it here in the form 



(43) 

with 

(44) 



(aq j ;q) j 



T = Y 



E T 



Hj{q) 



-2q 



OO 

u=l 



^+2 



E 

3=0 



(gq 3 \q)i 
(q~ 3 ;q)j 



H j(q)> 



v v qU+ 1 )' / 
~2q~ 



We will now study the functions Hj and propose to show that those of greater index contribute less significant 
terms in the asymptotic expansion of PA{q) near q — 1/2. In this way, a complete asymptotic expansion of 
the function PA(q), hence of its coefficients PA n , can be obtained. 

The main technique used here is that of Mellin transforms: we refer the reader to [T7] for details of the 
method. The principles are recalled in ^3.2.1 below. We then proceed to analyse the double sum T of (43) 
when q is real and q tends to 1/2. The corresponding expansion is fairly explicit and it is obtained at a 
comparatively low computational cost in §3.2.2 We finally show in §3.2.3| that the expansion extends to a 
sector of the complex plane around q — 1/2. 
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3.2.1. Principles of the Mellin analysis. Let f(x) be a complex function of the real argument x. Its 
Mellin transform, denoted by f*(s) or .M[/], is defined as the integral 

/•OO 

(45) M [f](s) = f*(s) := / f(x)x s - 1 dx, 



Jo 

where s may be complex. It is assumed that f(x) is locally integrable. It is then well known that if / satisfies 
the two asymptotic conditions 

f(x) = 0(x a ), f{x) = Otf), 

x-^0 x— f+oo 

with a > f3, then /* is an analytic function of s in the strip of the complex plane, 

-a < < -j3, 

also known as a fundamental strip. Then, with c any real number of the interval (— a, — /3), the following 
inversion formula holds (see [38, §VI.9] for detailed statements): 

i />c-f-zoo 

(46) fi x )= f*is)x~ s ds. 

ZlTT J c—ioQ 

There are then two essential properties of Mellin transfoms. 

(Mi) Harmonic sum property. If the pairs (A, /it) range over a denumerable subset of R x R >0 then one 
has the equality 



(47) M 



(a,aO 



/*(*) 




That is to say, the harmonic sum ^ A/(/xx) has a Mellin transform that decomposes as a product 
involving the transform of the base function (/*) and the generalized Dirichlet series (X^A/i ~ s ) 
associated with the "amplitudes" A and the "frequencies" fi. Detailed validity conditions, spelled 
out in [17] . are that the exchange of summation (^, in the definition of the harmonic sum) and 
integral ( f, in the definition of the Mellin transform) be permissible. 
(M2) Mapping properties. Poles of transforms are in correspondence with asymptotic expansions of the 
original function. More precisely, if the Mellin transform F* of a function F admits a meromorphic 
extension beyond the fundamental strip, with a pole of some order m at some point sq € C, with 
3?(so) < then it contributes an asymptotic term of the form P(logx)x _s ° in the expansion 
of Fix) as x — » 0, where P is a computable polynomial of degree m — 1. Schematically: 

(48) F*is) s4so =► Fix) _ o POog*)*— = Res (rW*" - )^ 

Detailed validity conditions, again spelled out in [T7], are a suitable decay of the transform F*(s), 



as 5(s) — >• ±00, so as to permit an estimate of the inverse Mellin integral (46) by residues - in (48 1, 
the expression is then none other than the residue of /*(s)x~ s at s — sq. 

The power of the Mellin transform for the asymptotic analysis of sums devolves from the application of 
the mapping property (M 2 ) to functions Fix) = ^Xfipx) that are harmonic sums in the sense of (Mi). 



Indeed, the factorization property (47) of (Mi) makes it possible to analyse separately the singularities 
that arise from the base function (via /*) and from the amplitude-frequency pairs (via ^A/i~ s ); hence an 
asymptotic analysis results, thanks to (M2). 



3.2.2. Analysis for real values of q — > 1/2. Our purpose now is to analyse the quantity T of (43) with 



q < 1/2, when q — > 1/2. This basically reduces to analysing the quantities Hjiq) of (44). Our approach 
consists of setting t = 1 — 2q and decoupling the quantities t and q. Accordingly, we define the function 



(49) h j it) = h j it;q,v):=q 



2^ NT 



An instance of such a decoupling technique appears for instance in Dc Bruijn's reference text [6] p. 27]. 
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so that 



with the definition (44 1. We shall let t range over K>o but restrict the parameter q to a small interval 
(1/2 — eo, 1/2 + eo) of K and the parameter v to a small interval of the form (3/2 — ei,3/2 + eil), since 
v(l/2) — 3/2. We shall write such a restriction as 

1 3 

with the understanding that eo, t\ can be taken suitably small, as the need arises. Thus, for the time being, 
we ignore the relations that exist between t and the pair q, v, and we shall consider them as independent 
quantities. 

As a preamble to the Mellin analysis, we state an elementary lemma 

Lemma 6. Let q be restricted to a sufficiently small interval containing 1/2 and v to a sufficiently small 
interval containing 3/2. Each function hj(t) defined by (49) satisfies the estimate 



(50) 



hAt) = 0[- 



0(1) ifj>l 

0(t~ 7 ) if j = 0, with 7 = 



\ogv 

log(iAz) ' 



For 7, we can also adopt any fixed value larger than log 2 (4/3) = 0.415, provided q and v are taken close 
enough to 1/2 and 3/2, respectively. 

Proof. Behaviour as t — > +oo. The inequality (l + tq^ l/ ^ 2 )^ 1 < t^ 1 q v+2 implies by summation the inequality 

hj(t) < q-H- 1 V v q> v q v = O (~ Y t -> +0O, 
i/=l ^ ' 

given the convergence of the geometric series for v « 3/2 and q ~ 1/2. 

Behaviour as t — > 0. First, for the easy case j > 1, the trivial inequality (1 + tq^^^ 2 )^ 1 < 1 implies 



hj(t) = o(^2(vq^A=0(l), *->0. 



Next, for j = 0, define the function 



i>b(t) := -2 + 



log(iA) 



log(l/g)' 

so that tq- y - 2 < 1, if v < v (t), and tq-"- 2 > 1, if 1/ > i/ Q (t). Write = Er/ + E„>„ • The sum 
corresponding to ^ > is bounded from above as in the case of t — > +oo, 



E 1 + v tq -,-2 ^ E «"t-v +2 = o(*- i («gr) = o(r i (« ? r), *->o, 



u>v (t) 



u>v (t) 



and the last quantity is 0(i 7 ) for 7 = (logw)/log(l/g). The sum corresponding to v < v$ is dominated by 
its later terms and is accordingly found to be 0(t~ 7 ). The estimate of ho(t), as t — > 0, results. □ 

We can now proceed with a precise asymptotic analysis of the functions hj(t), as t — » 0. Lemma [6] implies 
that each hj(t) has its Mellin transform h*(s) that exists in a non-empty fundamental strip left of 5J(s) = 1. 
In that strip, the Mellin transform is 



M[h s (t)] = q- 2 M 

( 51 ) = q- 2 M 

= q 



1 



l + t 
1 



E( u ? J ) I/ (9 " 2 ) S ) (by the Harmonic Sum Property (Mi)) 

1/= 



\V=1 



l+t 

2 ^ v <i 



1 - 

j+3s 



sin 7ts 1 — vqi +s 



(by summation of a geometric progression) 
(by the classical form of M[(l + t)^ 1 ]). 
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The Mellin transform of (1 + t), which equals 7r/sin(7rs), admits < 5i(s) < 1 as the fundamental strip, so 
this condition is necessary for the validity of (51 1. In addition, the summability of the Dirichlet series, here 
plainly a geometric series, requires the condition |tiq : ' +s | < 1; that is, 

log v 



> -3 + 



log 1/q' 



In summary, the validity of (51 ) is ensured for s satisfying 



A < Sft(s) < 1, with A := max ( 0, -j + U) '"' 



log 1/q 



Lemma 7. For q w 1/2 and v s» 3/2 restricted as in Lemma^ the function hj(t) admits an exact 
representation, valid for any t € (0,g -3 ), 



+ .?"-7nYl„„ *•> , „-2T/ i\r v 1 t r 



(52) ^(t) = (-iy ] p-m(io gl/q t) + q - 2 J2(-^ r 7 ; . 

log 1/q lq ^ l-f? 3 ' 

-ffere, 



/ \ logw 

7 = 7(9) — 



logl/g 

so that 7 rs log 2 § = 0.415, iu/jen q rs |; f/ie quantity H{u) zs an absolutely convergent Fourier series, 



(53) n(«):=^ ft e- a ", 
with coefficients pk given explicitly by 

7T 

(5 ' Pk = sin(7r7 + 2i^ 2 /(logl/9)) ' 

Observe that the pk decrease geometrically with k. For instance, at q = 1/2, one has 

(55) Pk = O ( e - 2fc - 2 /io g2 ^ = Q ( 4 28 . 1Q -i3j* i 

as is apparent from the exponential form of the sine function. Consequently, even the very first coefficients 
are small: at q — 1/2, typically, 

= |p_i| = 2.69 • 1(T 12 , \ P2 \ = \ P -2\ = 1.15 • 1(T 24 , |p 3 | = |p_ 3 | = 4.95 • 10~ 37 . 

Proof. We first perform an asymptotic analysis of hj (t) as t — > + . This requires the determination of poles 
to the left of the fundamental strip of h*(s), and these arise from two sources. 

- The relevant poles of ir/ sin ns are at s = 0, — 1, —2, . . .; they are simple and the residue at s = — r is 

(-iy. 

— The quantity (1 — vqi +s ) 1 has a simple pole at the real point 

56) o-o ■■= -J + ; — 

log 1/q 

as well as complex poles of real part o~o, due to the complex periodicity of the exponential function 
(e t+2 " = e*). The set of all poles of (1 - vq^ 8 )' 1 is then 

2ikir 

°~0 + 1 Tl~i G ^ 

log 1/q 



The proof of an asymptotic representation (that is, of (52 1, with '~' replacing the equality sign there) is 
classically obtained by integrating h*(s)t~ s along a long rectangle with corners at —d — iT and c + iT, where 
c lies within the fundamental strip (in particular, between and 1) and d will be taken to be of the form 
— m — i, with m £ Z> , and smaller than — j + 7. In the case considered here, there are regularly spaced 
poles along 3t(s) = — j + 7, so that one should take values of T that are such that the line 9(s) = T passes 
half-way between poles. This, given the fast decay of 7r/sin7rs as increases and the boundedness of 

the Dirichlet series (1 — vq : > +s )~ 1 along 5(s) = ±T, allows us to let T tend to infinity. By the Residue 
Theorem applied to the inverse Mellin integral (46), we collect in this way the contribution of all the poles 
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at — j + 7 + 2ifc7r/(log 1/q), with k £ Z, as well as the to + 1 initial terms of the sum in (52). The resulting 
expansion is of type (52) with the sum ~J2 r truncated to to + 1 terms and an error term that is 0{t m+1 / 2 ). 

In general, what the Mellin transform method gives is an asymptotic rather than exact representation 
of this type. Here, we have more. We can finally let to tend to infinity and verify that the inverse Mellin 
integral (46) taken along the vertical line 3t(s) = —to — | remains uniformly bounded in modulus by a 
quantity of the form ct m q~ 3m , for some c > 0. In the limit m 
tq~ 3 < 1), and the exact representation (52) is obtained. 



+oo, the integral vanishes (as long as 

□ 



We can now combine the identity provided by Lemma [7] with the decomposition of the generating 
function PA(q) as allowed by Equations (43) and (44), which flow from Proposition [3] We recall that 
Hj(q) = hj(t;q,v(q)). 

Proposition 4. The generating function PA(q) of prudent polygons satisfies, for q in a small enough 
interva^ of the form (1/2 — e, 1/2) (for some e > 0), the identity 

(a;q) oa (v;q) c 



(57) 



PA{q)=D{q)-q'A(q)- 



-T(q), 



(q;q) aD (av;q) 

C 

where the notations are those of Proposition^ and the function T(q) admits the exact representation 



(58) 



T(q) = (1 - 2 ff )-T • n 



log(l - 2g) 
logl/g 



U(q) + V(q), 1 



logw 
logl/g' 



withH(u) given by Lemma^ Equations (53) and (54). Set 

t = \-1q. 

The "singular series" U(q) is 



(59) U(q) = 

and the "regular series" V(q) is 

(q;q)oc q~'' 



vq 



37-2 (_ 



(-? _1 *;9)c 



(60) 



V(q) = 



log 1/q (-aq H^q)^ ' 
-2(9; 9)oo (av;q) c 



(a;<?)oo l + q 2 t 



(a;q) oc (v;q) c 



logu 
logl/<?' 



00 1 -1 -1 \ 

fl V 9 <l)r f _2.\ 

/ — i \ — (-aq t) 



Proof. We start from T{q) as defined by (J43J) . The g-binomial theorem is the identity [2TJ §1.3] 

(6z;q) c 



(61) 



(z;q)c 



Now consider the first term in the expansion ( 52 ) of Lemma [7] Sum the corresponding contributions for all 

(efficient 

,3 7 -2^ [aq-i;q)i 



values of j > 0, after multiplication by the coefficient ^q-j.'qy 1 m accordance with (43) and (44). This gives 



U(q) = 



vq ■ \- 

hgl/q^\q-3;q) 3 



i-qhy 



A simple transformation of type ( |42| ) finally yields 

(aq- j ]q)j. 2 . ■ _ ^ (a -1 ?; q)j 



3=0 vy j=0 



(-a<T 2 ^ 



which provides the expression for U{q) of the singular series, via the q-binomial theorem (61) taken with 
z = —at and 6 = a~ 1 q. 

Summing over j in the second term in the identity ( |52[ ) of Lemma [7j we have 



n 9 ) = <r 2 £Hr 2 i)T 



r=0 



(eg J ;q) 3 vq> r 
(l~ j \Q)j 1-vqi- 



Numerical experiments suggest that in fact the formula ( |58| remains valid for all q £ (0, 1/2). 
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Now, the Mittag-Leffler expansion ( 38 ) associated with Lemma [5] can be put in the form 
An application of this identity to V{q), with z = vq~ r , shows that 

v( q ) = ii-q-^ty ( _ x ) , 

( a ; 9)00 V N ;?)oo / 

which is equivalent to the stated form of V(q). Note that this last form is a g-hypergeometric function of 
type 2^i ; see [21]. 

So far, we have proceeded formally and left aside considerations of convergence. It can be easily veri- 
fied that all the sums, single or double, involved in the calculations above are absolutely (and uniformly) 
convergent, provided t is taken small enough (i.e., q is sufficienty close to 1/2), given that all the involved 
parameters, such as a,u,v, then stay in suitably bounded intervals of the real line. □ 

3.2.3. Analysis for complex values of q — » 1/2. We now propose to show that the "transcendental" 
expression of PA(q) provided by Proposition [4] is actually valid in certain regions of the complex plane that 
extend beyond an interval of the real line. The regions to be considered are dictated by the requirements of 
the singularity analysis method to be deployed in the next subsection. 

Definition 3. Let 8q be a number in the interval (0, ir/2), called the angle, and r$ a number in M>o, called 
the radius. A sector (anchored at 1/2) is comprised of the set of all complex numbers z = 1/2 + re 1 ® such 
that 

< r < r and 9 < 6 < 2tt - 9 . 

We stress the fact that the angle should be strictly smaller than tt/2, so that a sector in the sense of the 
definition always includes a part of the line 3?(s) = 1/2. The smallness of a sector will be measured by the 
smallness of tq. That is to say: 

Proposition 5. There exists a sector Sq (anchored at 1/2), of anglinOo < tt/2 and radius ro > 0, such that 



the identity expressed by Equations (57 1 and (58) holds for all q £ Sq 



Proof. The proof is a simple consequence of analytic continuation. We first observe that an infinite product 
such as (c; q)^ is an analytic function of both c and q, for arbitrary c and \q\ < 1. Similarly, the inverse 
l/(c; 9)00 is analytic provided cq> 7^ 1, for all c. For instance, taking c = a where a = a(q) = q 2 / (l—q+q 2 ) and 
noting that a(l/2) = 1/3, we see that l/(a; q) is an analytic function of q in a small complex neighbourhood 
of q = 1/2. This reasoning can be applied to the various Pochhammer symbols that appear in the definition 
of T(q), U{q), V{q). Similarly, the hypergeometric sum that appears in the regular series V(q) is seen to be 
analytic in the three quantities a f=a 1/3, v w 3/2, and t = 1 — 2q w 0. In particular, the functions U(q) and 
V{q) are analytic in a complex neighbourhood of q = 1/2. 
Next, consider the quantity 

(l-2g)-T = exp(-7log(l-2?)). 
The function 7 = -f(q) is analytic in a neighbourhood of q = 1/2, since it equals (logu)/(logl/g). The 
logarithm, log(l — 2q), is analytic in any sector anchored at 1/2. By composition, there results that (1 — 2g)~ 7 



is analytic in a small sector anchored at 1/2. It only remains to consider the II factor in (57). A single 

Iq, is also a 

Iog(l-2g)' 



Fourier element, p^e 2lk7TU , with u — \og x i q t and t = 1 — 2q, is also analytic in a small sector (anchored 
at 1/2), as can be seen from the expression 



(62) p k e-^™=p k exp [-2ikir 

V lo s 1 /q 

Note that, although 3t(log(l — 2q)) —> 00 as q —> 1/2, the complex exponential exp(2i/c7rlog 2 (l — 2q)) remains 
uniformly bounded, since 9(log(l — 2q)) is bounded for q in a sector. Then, given the fast geometric decay 
of the coefficients p k at q = 1/2 (namely, p k = 0(e~ 2kjr / lo e 2 ); cf. (54)), it follows that II(log 2 t) is also 



7 A careful examination of the proof of Proposition |5| shows that any angle Qq > 0, however small, is suitable, but only the 
existence of some 6>o < tt/2 is needed for singularity analysis. 
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analytic in a sector. A crude adjustment of this argument (see (71 1 and (72) below for related expansions) 
suffices to verify that the geometric decay of the terms composing (62) persists in a sector anchored at 1/2, 
so that II(log 1 / g t) is also analytic in such a sector. 

Finally, the auxiliary quantities D(q), A(q) are meromorphic at q = 1/2, with at most a double pole there; 
in particular, they are analytic in a small enough sector anchored at 1/2. We can then choose for iSo a 
small sector that satisfies this as well as all the previous analyticity constraints. Then, by unicity of analytic 



continuation, the expression on the right-hand side of (57), with T(q) as given by (58), must coincide with 
(the analytic continuation of) PA(q) in the sector So. □ 

3.3. Singularity analysis and transfer. If we drastically reduce all the non-singular quantities that occur 
in the main form (57) of Proposition [4] by letting q — » 1/2, we are led to infer that PA(q) satisfies, in a sector 
around q — 1/2, an estimate of the form 



(63) 

where 

(64) 



PA(q) = £ (1 - 2q)-^ 2 n(log 2 (l - 2q)) + O ((1 - 2 9 )" 



3/2 



Co 



1 (1/3; 1/2)00(3/2; l/2) c 

16 1 1 ' (1/2; 1/2)^(1/2; l/2) c 



U(l/2) = 



70 :=log 2 (3/2), 
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91og2' 



and U{q) is the singular series of (59). Let us ignore for the moment the oscillating terms and simplify H(u) 



to its constant term po, with pk given by (54). This provides a numerical approximation PA(q) of PA(q). 
With the general asymptotic approximation (derived from Stirling's formula) 

1 



(65) 



[«"](! -2 9 ) 



— A 



,A-1 



>+oo r(A)' 



A 



'<0) 



it is easily seen that \q n ]PA(q) is asymptotic to the quantity Ko2 n n 9 of Equation (18) in Theorem [2] which 
is indeed the "principal" asymptotic term of PA n = [q n ]PA(q) 7 where g = 70 + 1 = log 2 3. 

A rigorous justification and a complete analysis depend on the general singularity analysis theory |201 
Ch. VI] applied to the expansion of PA(q) near q = 1/2. We recall that a /^.-domain with base 1 is defined 
to be the intersection of a disc of radius strictly larger than 1 and of the complement of a sector of the form 
— 9q < arg(z — 1) < 6*o for some #0 € (0, 7r/2). A A-domain with base p is obtained from a A-domain with 
base 1 by means of the homothetic transformation z 1— >• pz. Singularity analysis theory is then based on two 
types of results. 

(Si 

p. 381]. In the case of the scale (1 



Coefficients of functions in a basic asymptotic scale have known asymptotic expansions [20jj Th. VI. 1, 

A , the expansion, which extends (65), is of the form 



m-z) 



— A 



,A-1 



n— ^+00 



1 + 



\ - e fc 



AeC\Z> , 



fc>i 



where e/- is a computable polynomial in A of degree 2k. Observe that this expansion is valid for 
complex values of the exponent A, and if A = a + ir, then 



,A-1 



-lgirlogn 



Thus, the real part (a) of the singular exponent drives the asymptotic regime; the imaginary part, 
as soon as it is nonzero, induces periodic oscillations in the scale of log n. A noteworthy feature is 
that smaller functions at the singularity z = 1 have asymptotically smaller coefficients. 
(S2) An approximation of a function near its singularity can be transferred to an approximation of coef- 
ficients according to the rule 

f(z) = 0((l-z)- A ) [z n ]f(z) = O^- 1 ). 

The condition is that f(z) be analytic in a A-domain and that the O-approximation holds in such a 
A-domain, as z 1; see Th. VI. 3, p. 390]. Once more, smaller error terms are associated with 
smaller coefficients. 



Equipped with these principles, it is possible to obtain a complete asymptotic expansion of [q n ]PA(q) once 
a complete expansion of PA(q) in the vicinity of q = 1/2 has been obtained (set q = z/2, so that z « 1 
corresponds to q = 1/2). In this context, Proposition [3] precisely grants us the analytic continuation of 
PA(q) in a A-domain anchored at 1/2, with any opening angle arbitrarily small; Proposition [4j together 
with Proposition [5] describe in a precise manner the asymptotic form of PA(q) as q — > 1/2 in a A-domain 
and it is a formal exercise to transform them into a standard asymptotic expansion, in the form required by 
singularity analysis theory. 

Proposition 6. As q — > 1/2 in a A-domain, the function PA(q) satisfies the expansion 

(66) PA(q) ~ ^^ + i ? ( (7)+ ^ ( i_2 g )-^- 2 +^(log(l-2g)) £ n(^(log 2 (l-2 (? )), 7o = log 2 §. 
Here R(q) is analytic at q = 1/2 and eac/i 11^'^ (m) is a Fourier series 

kez 

with a computable sequence of coefficients p^'^ ■ 
Proof. From Proposition |4j we have 

(67) PA(q) = PA ICS (q) + PA sins (q), 

where the two terms correspond, respectively, to the "regular" part (involving the regular series V(q)) and 
the "singular part" (involving the singular series U(q) as well as the factor (1 — 2q)~ 7 and the oscillating 
series) . 

Regarding the regular part, we have, with the notations of Proposition |4j 

(68) PA^(q) = D(q) - q *A(q) { °[ q \~ { ": q) °° V(q). 

We already know that A(q) and D(q) are meromorphic at q = 2 with a double pole, while V(q) and the 
Pochhammer symbols are analytic at q = 1/2. Thus, this regular part has at most a double pole at q = 1/2. 
A simple computation shows that the coefficient of (1 — 2q)~ 2 reduces algebraically trivially - in the sense 
that no g-identity is involved - to 0. Thus, the regular part involves only a simple pole at q = 1/2, as is 



expressed by the first two terms of the expansion (67), where R(q) is analytic at q = 1/2. (Note that the 



coefficient of (1 — 2q)~ 1 is exactly 1, again for trivial reasons.) 

The singular part is more interesting and it can be analysed by the method suggested at the beginning 
of this subsection. Whenever convenient, we freely use the abbreviation t = 1 — 2q. The function 7(g) = 
(logi>)/(logl/g) is analytic at q = 1/2, where 

, x 3 log 3 , i . 

7(9) = l0 S 2 2 +2 (loi2F ( ^^ + --- 

= 0.58496 + 4.5732(q - §) + 16.317(q - \f + 39.982(q - i) 3 + 86.991(q - \ f + • • • . 

The function (1 — 2q) -7 can then be expanded as 

{l-2q)-~< {q) = {l-2q)-^e- ( -^-^ lo& \ with 7o = 7 (l/2) = log 2 | 



= (l-2q)-^' [l + —^—t\ogt + t 2 P 2 (\ogt)+t 3 P 3 (\ogt) 



i69) „ „ W.n (, , log 3 

v (log2) 2 ' 

for a computable family of polynomials P%, P3, . . ., where degP^ = I and Pi(0) — For instance, we have, 
with y := \ogt: 

(1 - 2q)- (7 ~ 7o)logt = 1 + 2.28% + i 2 (-4.07y + 2.61y 2 ) + t 3 (4.99y - 9.32y 2 + 1.99?/ 3 ) + • • ■ . 



The singular series U(q) of (59) is analytic at q = 1/2 and its coefficients can be determined, both numerically 
and, in principle, symbolically in terms of Pochhammer symbols and their logarithmic derivatives (which lead 
to g-analogues of harmonic numbers). Numerically, they can be estimated to high precision, by bounding the 
infinite sum and products to a finite but large value. (The validity of the process can be checked empirically 



by increasing the values of this threshold, the justification being that all involved sums and products converge 
geometrically fast - we found that replacing +00 by 100 in numerical computations provides estimates that 
are at least correct to 25 decimal digits.) In this way, we obtain, for instance, the expansion of the function 
V(q), which is of the form (t = 1 — 2q) 



9.97t + 21.5i 2 + 35.8t 3 + 51.9t 4 



Finally, regarding II(u) taken at u = log 1 / q (l — 2q), we note that the coefficients pk of (54 1 can be 
expanded around q = 1/2 and pose no difficulty, while the quantities e 2lk7TU can be expanded by a process 
analogous to (69). Indeed, we have 

(71) 



Pk =Pk(q) 



exp (1 + e\{k)t + e 2 (k)t 2 - 
log 1 / ? (l — 2q), one has 



sin(7T7o + 2ikn 2 / log 2) 
where the e& only grow polynomially with k. Also, at u 

(72) e 2lk7TU = (1 - 2 g) 2 ^/ log2 exp [2zfc^log 2 t ( 9l t + g 2 t 2 + ■■■)], 

where the coefficients gj are those of (logl/g) -1 — (log2)~ 1 expanded at q = 1/2 and expressed in terms 
of t = 1 - 2q. 



We can now recapitulate the results of the discussion of the singular part: from (70 1, (71 1, and (72), we 
find that the terms appearing in the singular expansion of PA(q) are of the form, for j = 0, 1, 2, . . ., 

(1 - 2q)-^- 2 V (logt) e t 2lk */ los2 , 

with I such that < I < j and k <E Z. The terms at fixed j,£ add up to form the Fourier series H^ ,£ \ 
whose coefficients exhibit a fast decrease with |fc|, similar to that encountered in (55). Consequently, a finite 
version of (66) at any order holds, so that the statement results. □ 

With this last proposition, we can conclude the proof of Theorem [2] 

Proof (Theorem^). The analytic term R(q) in ( [66] ) leaves no trace in the asymptotic form of coefficients. 
Thus the global contribution of the regular part to coefficients PA n reduces to 2™ (with coefficient 1 and no 
power of n), corresponding to the term (1 — 2q)~ l in (66). 

The transfer to coefficients of each term of the singular part of (66) is permissible, given the principles of 
singularity analysis recalled above. Only an amended form allowing for logarithmic factors is needed, but this 
is covered by the general theory: for the translation of the coefficients of the basic scale (1 — z)~ A log fc (l — z), 
see [201 P- 387]. From a computational point of view, one may conveniently operate [201 Note VI. 7, p. 389] 
with 



[z"](l-z)- A (log(l-z)) fc = (-l)' 



d T{n + X) 



<9Ar(A)r(n + l)' 

then replace the Gamma factors of large argument by their complete Stirling expansion. 

We can now complete the proof of Theorem[2] It suffices simply to retain the terms corresponding to j = 
in (66), in which case the error term becomes of the form O ((1 — 2q) _7o_1 log(l — 2g)), which corresponds 
to a contribution that is 0(n 7 ° logn) = 0(n 9 ~ 1 logn) for PA n . 

Next, regarding the Fourier element of index k = 0, the function-to-coefficient correspondence yields 

-1 



(1 - 2q) 



-70-2 



n 
r (7o 



l + O 



Thus, the coefficient kq in ( pL9[ ) has value (cf (63) and (64)) given by 

1 



«o = Co • Pol 



9=1/2 ' r( 7o + 2)' 



7o 



log 3/2 
log 2 ' 



with Co as m (64). This, given the form (|54j) of pk at k = 0, is equivalent to the value of n stated in 
Theorem [2] (where g :— 70 



1 = log 2 3). 

For a Fourier element of index k £ Z, we have similarly 

,,70 + l+iXfc 

(1 - 2q)-^- 2 - lXk => —, , (1 

r( 7o 1 - 1 - 



I- iXk) 

23 



0(l/n)) 



where 



Xk 



2kn 
log2' 



FIGURE 6. Left: The difference (PA n — Q.§)2~ n n~ 9 against log 2 n, where fie is the six-term exten- 
sion of (74 1, for n = 60, . . . ,500; the plot reveals about three periods of the oscillating component 
in the asymptotic expansion of PA n . Right: The corresponding plot, for values up to n = 6000. 



The diagrams confirm the presence of oscillations that, asymptotically, have amplitude of the order 
of 10 -9 . 



We finally observe that 



n 7o+lgiXfc log" 



so that all the terms, for k 6 Z, are of the same asymptotic order (namely, O(n 7o+1 )) and their sum 
constitutes a Fourier series in logn. The Fourier coefficient then satisfies, from the discussion above: 

Thus finally, with g = 70 + f : 

(73) K = - ft (i-|2-j)(i-|2-j) 

1 ' k 91og2sin(7rg + 2ifc7r 2 /log2)r( 5 + l + 2iA:7r/log2)| = -l ) (1 - §2-J) 2 

This completes the proof of Theorem [2j □ 

The same method shows the existence of a complete asymptotic expansion for PA n . 
Theorem 3. The number of 3-sided prudent polygons satisfies a complete asymptotic expansion, 

PA n ~ 2" + 2" • n 9 (So,o + - Qogn ■ ^1,1 + 5 1>0 ) + \ (log 2 n • S 2 , 2 + logn • S a ,i + 3 2 , ) ■ 

where Sj^ is an absolutely convergent Fourier series in logn. 

The non-oscilating form obtained by retaining only the constant terms of each Fourier series is computed 
by a symbolic manipulation system such as Maple or Mathematica in a matter of seconds and is found 
to start as 



^5 
2" 



0.1083842947 ■ n 9 



+ (-0.3928066917L + 0.5442458535) • n 3 ' 1 
( 74 ) + (0.2627062704 1? + 0.6950193894 L + 0.6985601031) • n 9 ' 2 

+ (0.08310555463 L 3 - 0.02188678892 1? - 1.570478457 L - 1.18810811075202) • n 9 " 3 

+ (0.06722511293 L 4 + 0.05494834609 L A - 3.297513638 L 2 - 4.663711650 L - 4.156441653) ■ n 9 " 4 , 

where L = logn. In principle, all the coefficients have explicit forms in terms of the basic quantities 
that appear in Theorem [2] (augmented by derivatives of g-Pochhammer symbols at small rational values) . 



However, the corresponding formulae blow up exponentially, so that we only mention here the next coefficient 



0.39280 ... in (74), whose exact value turns out to be 



log 3 
log 2 

Figure [6] displays the difference between PA n and the six-term extension fig of ( 74 ) . It is piquant to note that 
all the terms given in Equation ( 74 ) are needed in order to succeed in bringing the fluctuations of Figure [6] 
out of the closet. As a matter of fact, before the analysis of Theorem [3] was completed, the authors had 
tried empirically to infer the likely shape of the asymptotic expansion of PA n , together with approximate 
values of coefficients, from a numerical analysis of series. The success was moderate and the conclusions 
rather "unstable" (leading on occasion to heated debates between coauthors). This state of affairs may 
well be present in other yet-to-be-analysed models of combinatorics. It is pleasant that, thanks to complex 
asymptotic methods, eventually . . . everything nicely fits in place. 

4. Conclusions 

Classification of prudent polygons. Our first conclusion is that the present study permits us to advance 
the classification of prudent walks and polygons: the generating functions and their coefficient asymptotics 
are now known in all cases up to 3-sided (walks by length; polygons by either perimeter or area). Functional 
equations are also known for 4-sided prudent walks and polygons, from which it is possible to distill plausible 
estimates. We can then summarize the present state of knowledge by the following table (compare with Q 
and ^). 



Polygons (area) Generating function Asymptotic number References 



2-sided rational 2™ this paper, [ 2.2 



3-sided non-holonomic C-${n) • 2 n n ° S2 this paper, Th 



4-sided functional equation C*4(n) • 2 n n 1+log2 3 (?) this paper, [ 2.4 



Beaton et al. p]. 



(The oscillating coefficient C^lri) is expressible in terms the Fourier series k(u) of Theorem[2j) The numerical 
data relative to 4-sided polygons enumerated by area (last line) are from pQ: in all likelihood, the dominant 
singularity remains at 1/2 and the critical exponent is 1 + log 2 3; that is, one more than the corresponding 
exponent for 3-prudent polygons. Examination of subdominant terms also suggests that the number of 
4-polygons, once divided by n, satisfies an expansion of the form obtained in Theorem [3] for the 3-prudent 
counterparts. The "mean" amplitude is probably about 0.033, and, under the circumstances, there is little 
doubt that minuscule oscillations (rendered by Ci(n)) must also be present. The foregoing analysis of the 

3- prudent case then at least has the merit of pointing towards the type of singularity to be expected for 

4- prudent polygons as well as, possibly, to methods of attack for this case. 

Methodology. The generating function of 3-prudent polygons has been found to be a g-hypergeometric 
function, with the argument and parameters subject to a rational substitution. The methods developed 
here should clearly be useful in a number of similar situations. Note that the asymptotic enumeration of 
prudent walks and polygons by length and perimeter is in a way easier since the dominant singularity is 
polar or algebraic. (Bousquet-Melou [3] however exhibits an interesting situation where the complete singular 
structure has a complex geometry.) 

Estimates involving periodic oscillations are not unheard of in combinatorial asymptotics [201 HH l2"9"] . 
What is especially interesting in the case of 3-sided prudent polygons is the pattern of singularities that 
accumulate geometrically fast to the dominant singularity. This situation is prototypically encountered in 
the already evoked problem of the longest run in strings: the classical treatment is via real analytic methods 
followed by a Mellin analysis of the expressions obtained; see [35] • In fact, the chain 

(76) Coefficient asymptotics Mellin transform + Singularity analysis 

is applicable for moment analyses. For instance, the analysis of the expected longest run of the letter a in a 
random binary sequence over the alphabet {a, b} leads to the generating function [221 Ex. V.4] 

k 

$( Z ) = (1 - Z) V 7-t, 

\ J \ I i-2z + z k+1 ' 



k>0 
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to which the chain ( 76 ) can be applied. We could indeed recycle for our purposes some of the corresponding 
analysis of (SSJUO]) when discussing the location of poles of the generating function PA(z), in Section [3j 

Another source of similar phenomena is the analysis of digital trees [26l [34] , when these are approached 
via ordinary generating functions (rather than the customary exponential or Poisson generating functions). 
Typically, in the simplest case of node-depth in a random digital tree, one encounters the generating function 

1 2~ fe 

*(*) = 



1-z ^ 1 - z(l- 2- fe )' 



where the geometric accumulation of poles towards 1 is transparent, so that the chain ( 76 ) can once more 
be applied [19]. 

We should finally mention that "critical" exponents similar to g — log 2 3 surface at several places in 
mathematics, usually accompanied by oscillation phenomena, but they do so for reasons essentially simpler 



than in our chain (76). For instance, in fractal geometry, the Hausdorff dimension of the triadic Cantor 
set is 1/g, see {14 , while that of the familiar Sierpihski gasket is g, so that g occurs as critical exponent in 
various related integer sequences |18j . The same exponent g — log 2 3 ~ 1.58 is otherwise known to most 
students of computer science, since it appears associated to the recurrent sequence f n = n+ 3/i n /2j , which 
serves to describe the complexity of Karatsuba multiplication |27] (where, recursively, the product of two 
double-precision numbers is reduced to three single-precision numbers). In such cases, the exponent g is 
eventually to be traced to the singularity (at s — g, precisely) of the Dirichlet series 



which is itself closely related to the Mellin transforms of our Eq. (51 ). See also the studies by Drmota and 
Szpankowski [10] , Dumas [12] , as well as [18] for elements of a general theory. 

For the various reasons evoked above, we believe that the complex asymptotic methods developed in the 
present paper are of a generality that goes somewhat beyond the mere case of 3-sided prudent polygons. 
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